The Generalized DMPK equation revisited: A systematic derivation 
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The Generalized Dorokov-Mello-Pereyra-Kumar (DMPK) equation has recently been used to ob- 
tain a family of very broad and highly asymmetric conductance distributions for three dimensional 
disordered conductors. However, there are two major criticisms of the derivation of the Generalized 
DMPK equation: (1) certain eigenvector correlations were neglected based on qualitative arguments 
that can not be valid for all disorder, and (2) the repulsion between two closely spaced eigenvalues 
were not rigorously governed by symmetry considerations. In this work we show that it is possible 
to address both criticisms by including the eigenvalue and eigenvector correlations in a systematic 
and controlled way. It turns out that the added correlations determine the evolution of the Jaco- 
bian, without affecting the evaluation of the conductance distributions. They also guarantee the 
symmetry requirements. In addition, we obtain an exact relationship between the eigenvectors and 
the Lyapunov exponents leading to a sum rule for the latter at all disorder. 

PACS numbers: 73.23.-b, 71.30., 72.10. -d 
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I. INTRODUCTION 

A highly non-trivial, non-Gaussian distribution P(g) 
of conductances g for disordered three dimensional (3d) 
conductors in the large disorder limit has recently been 
predicted in Refs. [HQ. The predictions include, among 
others, a strong deviation from the expected log-normal 
distribution in the deeply insulating limit in agreement 
with numerical results obtained from tight binding An- 
derson model^, variance that grows with disorder con- 
sistent with numerical simulations-, and an asymme- 
try that changes sign near the metal-insulator transi- 
tion as the disorder is decreased from the deeply insu- 
lating limit. The distribution was obtained using the so 
called Generalized DMPK equation 7 -. While the orig- 
inal Dorokov-Mello-Pereyra-Kumar (DMPK) equation 8 
has been shown to be valid for quasi one dimensional 
(Id) systems^ where there is no Anderson metal-insulator 
transition, the Generalized DMPK equation has been 
claimed to be valid beyond quasi Id in the large disorder 
regime. Since the full distribution P(g) in 3d at strong 
disorder is beyond the scope of the conventional field the- 
ory framewor k 10 ' 11 , and since a broad and highly asym- 
metric distribution can have interesting consequences for 
the Anderson transition, it is important to critically eval- 
uate the validity of the Generalized DMPK equation. 

For a given iV-channel conductor of length L z and cross 
section L 2 , the 2N x 2N transfer matrix M relates the 
wavefunction on the left of the sample to that on the 
right via: = M^l- For the case of time reversal and 
spin rotation symmetry (orthogonal ensemble), M can 
be parameterized by the N x N unitary matrices u and 
v, and N x N diagonal matrix A: 



M 



u 
u* 



VTTx Vx 
■fx Vi + x 



v 
v* 



(1.1) 



the eigenvectors since they turn out to be the eigenval- 
ues and eigenvectors of a certain combination of MM\ as 
will be shown later. The |u 0& | 2 and |v a &| 2 can be thought 
of as the fraction of current in channel a that enters into 
channel b from left to right and from right to left, respec- 
tively, and 1/(1 + Aft) is channel fr's contribution to the 
conductance. M is a random matrix, and its probability 
distribution can be written as 

p(M)dn(M) = p(X, u, v)dn(X)dn(u)dn(v). (1.2) 

The measure dfi(X) can be written in terms of the so- 
called Jacobian J (A): 

N 

dtx{X) = J(X)d N X; J{X) = l[\X a - X b f (1.3) 

where the orthogonal symmetry discussed here requires^ 
;8=1. For future reference we'll define related functions: 

p(X,v) = / dfi(u)p(X,u,v), 



(1.4) 



p(X) — j d)i(u)dfj,(v)p(X,u,v), 
P(X) = p(X)J(X). 



The last one is the probability distribution of the eigen- 
values A, (i — 1,2,- ■ -N). One can obtain an evolu- 
tion equation for p(X) in the following manner. We 
add a small conductor of length 5L Z described by a 
transfer matrix M$l z = M' to a conductor of length 
L z described by Mr, z = M". Then taking advantage 
of the fact that the transfer matrices are multiplicative 
(M Ljs+ sl z = M — M"M'), an integral equation for the 
evolution of p(M) with length L z can be written down: 



We will loosely refer to A and it, v as the eigenvalues and 



p Lm +6L.(M)= / dii{M')p SLz {M')p L „{M h 



(1.5) 
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Integrating both sides over dfi(u) and dfi(v), we can write 
the evolution equation for p(X) as 

Pl z +sl z (A) = J dfi(u)d^(v)(pL z {X" ,u" ,v"))m' (1.6) 

where the angular bracket (• • • )m* indicates an averag- 
ing over M' . When the sample width is much smaller 
than the localization length £ (L « £), and the length 
is much larger than the width (L z » L), the matri- 
ces u and v are isotropically distributed^. In that case, 
p(X, u, v) only depends on A. To convert the integral 
equation to a differential equation, p(X") is expanded to 
second order in a Taylor series about A" = A, and the av- 
erage over M' and integrals over dfi(u) and dfi(v) can be 
done. The resulting evolution equation for P(X) (DMPK 
equation) is 



dP(X) 
dt 



2 N d 

NTi^ax- a Xa{1 + Xa) 



a=l 

d 

dXa 



Sir 2 



X n — Xh 



P(X) (1.7) 



where t = L z /l, and I is the mean free path. One can 
obtain the conductance distribution from P(X) via the 
Landauer formula^ 



H?)« / {[d N xp(x)s [9-J2t 

J 1=1 \ i 



1 



+ Xi 



(1.8) 



The V{g) obtained in this way has been shown to be 
valid for quasi Id systems at all disorder, with surprising 
features near the metal-insulator crossover regime^—. 

The first attempt towards relaxing the isotropy ap- 
proximation was made in Ref. 18]. It was assumed 
that p(X,u, v) depended only on A, but u and v were 
not isotropically distributed however. The p(X") was ex- 
panded out to second order about A" = A. The average 
over M' and integral over d/j(u) could be done and the 
integral over dfj,(v) was parameterized via arbitrary con- 
stants fii and fi2- Requiring that the resulting equation 
conserve probability necessiated a renormalization of the 
Jacobian from J to J 7 so that P(X) was redefined as 



N 



pw = n i a » 

a=£b 



(1.9) 



It also linked [i\ and /i2 together leaving their ratio 7 = 
as the only true degree of freedom, resulting in a 
one-parameter generalization of the DMPK equation: 



dP{X) 
dt 



N 



N 



2 " 2 ii:^,i + y 

a— 1 



d 

dXa 



E 



X a — Xb 



P(X). (1.10) 



The parameter \xi just renormalizes the mean free path, 
while the presence of 7 evinces a disorder-dependent re- 
pulsion between the eigenvalues. It was demonstrated in 
Ref. [lH that 7 must be 1 in the weak disorder isotropic 
limit, and asymptotically approach in the large disorder 
limit. 

Relaxing the assumption about p{X, u,w)'s exclusive 
dependence on A entirely, a more formal derivation of 
the Generalized DMPK equation was given in Ref. 
Going back to Eq. (11.61) , a model for the average over M' 
known to reproduce the metal-insulator transition was 
usedi^, and it was found that the integral over d[i(u) 
could be done exactly. The resulting p{X" , v") was Taylor 
expanded about A" = A, but it was assumed that an ex- 
pansion of v" about v was unnecessary, especially in the 
insulating regime of interest because in that regime the 
eigenvectors settle into length-independent values. The 
required integrals over d/i(u) were formally done using a 
mean field approxition using the following definitions: 



K ab = (k ab ) Lz = ^(\v aa \ 2 v ba \ 2 ) Lz 



and also 



lab 



2Ka 

K a , 



(1.11) 



(1.12) 



Here the angular bracket denotes an ensemble average. 
Requiring that the resulting equation obey probability 
conservation necessitated a renormalization of the Jaco- 
bian again so that P(X) was re-defined to be: 



N 



P(X) = Y[\X a -X b \^p(X). 

a^b 



(1.13) 



One could then write the resulting Generalized DMPK 
equation (GDMPK) as 



8P(X) 
dt 



N 







^ K aa ——X a {l + Xa) 



lab 



9X a j-f X a — X b 

b^a 



P(X). (1.14) 



As shown in Ref. [2|, the matrix elements K ab can be 
explicitly evaluated numerically and contain information 
not only about dimensionality but also the critical point 
in 3d. The quasi Id DMPK equation is recovered when 
lab — > 1. In the above derivation the eigenvector corre- 
lations were neglected under the assumption that they 
become small at very large disorder, restricting the va- 
lidity of Eq. (|1.14j) to the deeply insulating limit. Within 
this limit, the conductance distributions obtained ana- 
lytically from the solutions of Eq. (11.141) agree well with 
numerical results obtained from tight binding Anderson 
modet^. 

Despite the successes, there are a few major criticism of 
the derivation of Eq. (I1.14j) . The first one is that it is not 
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clear if and/or under what conditions the eigenvector cor- 
relations can be neglected. Second, it has heretofore been 
necessary to impose probability conservation as a sepa- 
rate constraint, which has required a 'renormalization' of 
J(A), and redefinition of P(A) seemingly at odds with Eq. 
(|1.4|) . Last, there is the apparent conflict with the known 
constraint that in the limit | A Q — A& | — > 0, the level repul- 
sion should be characterized by the symmetry parameter 
/3 = 1 as determined by the exponent of the Jacobian 
and not by the matrix 7 a ^, where 7 a fc <C 1 in the deeply 
insulating limit. In this work we address all of these 
criticisms. To address the first and second criticisms 
we include both eigenvalue and eigenvector correlations 
systematically in a controlled way up to order SL z /l. 
We will find that including the additional eigenvector 
correlations results ultimately in the same GDMPK in 
Eq. (|1.14p . but also restores the Jacobian to its origi- 
nal form in Eq. (jl.3l) , resulting in an equation which au- 
tomatically conserves probability. Thus the results for 
V(g) obtained earlier in Refs. [IQ turn out to be valid 
quite generally, beyond the deeply insulating limit. We 
address the last criticism by developing a self-consistent 
evolution equation for K a b and showing that j a i, — > 1 
as |A a — Afc| — > 0, satisfying the symmetry requirement. 
In addition, we obtain an exact relationship between the 
eigenvectors and the Lyapunov exponents v n , defined by 
the relation (A„) = exp[(2L z /l)u n ], leading to a sum rule 
for the Lyapunov exponents. 



II. THE GENERALIZED DMPK EQUATION 

We will start with Eq. (jl.5|) . and from there basically 
follow Ref . Q , except that we will not neglect the eigen- 
vector correlations. For the sake of completeness and 
comparison, we will rederive some of the results already 
obtained in Ref. 0. As in Ref. Q, we will go beyond the 
'isotropy' approximation of the original DMPK equation 
by keeping the u and v dependence of pl z in addition 
to the A dependence. Using Eq. (jl.ip . and equating M" 
with MM'- 1 , it then follows that u" =u- it" (A, A', vv' r ). 
The invariance of the measure then allows one£ to inte- 
grate over u, leading to an integral equation purely in 
terms of p(X, v), 



Pl z +sl z (X,v) = J dfi(u')dfj,{X')dfi(v')psL g (X' ,u' ,v') 
x PLz (X\v"). (2.1) 

Since pl z +sl z {X,v) is independent of u, it must be the 
case that ps L (u',v',X') is independent of u'. So we set 
v! = v'* for later convenience, and write 



PL,+5L m (X,v) = / d^(X')d^(v')ps Lz {^ '>"') 

x Pi ,(A",v"). (2.2) 

Next we write A" = A + SX" and v" = v + Sv" and 
Taylor expand p(X",v"). Thereafter we integrate both 



sides over dfi(v) and use integration by parts to write: 

Pl z +slA x ) =Pl z W 
d f 

+ E^ / d V( v )PLz&,v){5X"}M' 

1 d 2 f 

+ 2 E q\ d X / d ^ PL ^ v ^ 5X 'a 6X b)M> 

ab a 



9K h 

dv ab 



M 1 



E / d ^ V )PLA\ 

2^ / Mv)pl z A, v) ( a 



ab.cd 



dv ab dv cd 



M' 



(2.3) 



We now assume that the impurity averaged correla- 
tions < • • • > reach a limiting distribution that has a 
well-defined peak, possibly A - dependent, roughly coin- 
cident with their mean. For example in the weak disorder 
limit, the correlations have a peak around the quasi- Id 
distribution. As checked numerically in Ref. 0, this is 
always true provided the length L z is sufficiently large 
(L z 3> I) (and the disorder is not so large that £ ~ I ~ a 
where a is the lattice spacing) . Note that we do not need 
to know what the peaked distribution actually looks like, 
it is enough to assume that such a distribution exists. In 
that case, equation (|2.3[) can be rewritten as 



dL z 



where we define, 



/(*WA) + I>(A)^i 



(2.4) 



/(A) = 

9a (A) = 
hab(X) 



I 

I 

= SL Z 



d 2 SvX d 



2 i^,, \ 9v ab dv cd 

ab.cd 



E 

ab 



9Kt 

dv ab 



I 



2SL 



-(SK6K 



be ^ 



dvbc 



(2.5) 



Here the angular brackets (• • • ) represent average over 
both M 1 and v. Note that in the limit 5L Z — > 0, we only 
need to keep terms in the average that are at most linear 
in SL z /l. In Ref. 0, only the terms corresponding to the 
changes in A" were assumed to be important, neglecting 
the changes in 5v" . Since these also contain terms linear 
in 8L z /l, a systematic expansion needs to keep all five 
terms. All other terms in the series expansion beyond 
the second derivative (represented by (• • • ) in Eq. (|2.4|l ) 
will be higher order in SL z /l and therefore need not be 
considered. In order to evaluate the changes SX" and Sv", 
we would like a matrix whose eigenvectors and eigenval- 
ues are given by v and A respectively. Such a matrix, X, 



4 



is given below^S: 

X = -[M t M+(M t M) _1 -21] 

A little algebra, using Eq. is sufficient to demon- 

strate the last line. From here we form X" out of 
M" = MM'' 1 to calculate the perturbative corrections 
to A" and v". After some manipulations one obtains 



/'tA'V =v ] \v + v ] Wv 



(2.7) 



where the matrix W is given by (using vl = v'* and 
expanding out to linear order in A') 

W = w'U'wV + ~m/ t AVv t A + -Aw'+A'uV 
y/\(l + \)v*v' T VXv'vt 



- vv"'V>Jv'*v T v / X{l + A) 
+ vv^Vx 1 v'*v T Xv*v' T Vx 1 v'v^. 



(2. 



The perturbative changes due to M' are then given by 

\W m \ 2 



5X" 



E 



"t _ \ ^ 



SvJ„ = 



w„ 



A, 



A; 



1 



A„, — Xi 



iV 

1 rt 



\ " \ ~> Wjj vv Jn f 
Z—' Z—' X X X X " 



2 E ' ""' 



E 



Win t 



An ~~ Aj A n — Aj 



(2.9) 



Note that the eigenvector corrections so written are nor- 
malized to 2nd order. To make further progress, we follow 
Ref. [7| and use a very general model for the average over 
M'. Since A' — > as 5L Z — > 0, we use A' oc Sl^/I and the 
folowing 'building blocks': 



Kb) = (V'; b ) = KM = (v'aWcd) = 



N 



l\l I* I \ "Lz r 



^{\jKK V a* V aP v b5 v Wl M ' = -T- S <*l!> S orrSpS- (2-10) 



ab 



These are the building blocks used in the original DMPK 
equation and those suggested in Ref. [19( to describe the 
metal-insulator transition. This allows us to build a sys- 
tematic and controlled expansion of p^ z in powers of the 
small parameter SL z /l. All of the correlations in Eq. (I2.5[) 
can be written in terms of the following two results: 



8L, 



(1 + x n )s v 



5L Z s—^ 2 

—j- X a v mj3 v n p \v a p\ 



a/3 



(WijW mn ) M > 



SL, 



A in ^' 



i * V jm / f V ma V na V ia v ja 



where we make the definition 



Ay = ^Xiil + X^Xjil + Xj). 



(2-11) 



As is also evident in Eq. (12.51) . we need to take deriva- 
tives with respect to v a b and we use the following two 
basic equations for this, 



dv 



ab 



dvi 



(2.12) 



The latter is obtained by implicitly differentiating the 
relationship vnv — 1. Finally, we need to then calculate 
averages over v. To facilitate this, we make a random 
phase assumption, that averages over products of u's will 
be zero unless their phases cancel. With these in hand, 
and after lengthy algebra we find that /, g, and h of 
Eq. (|2.5p can all be written in terms of the K matrix, 



/(A) = £ [A'- 



k^m 

+ E 



Am — Afe 

A m (l + A m ) 



2K n 



(A m — Afc)(A m — X n ) 



[Kn 



Ja(X) = {l + 2X a )K aa + Y, 



2A Q (1 + A a ) 
A a — Afe 



2K n 



[K aa — K a k] 



hab(X) = A a (l + X a )K aa 5 a b- 



(2.13) 



Putting these results into Eq. (|2.4|) and using the defi- 
nition of P(X) from Eq. (|1.4p we find that we can write 
Eq. (|2.4I) in full generality as 



dP(X) 
dt 



N 







= ^X aa — A (1 + A Q ) 



dX a 

dX a {-f X a - X b 

b^a 



lab 



P(A). (2.14) 



As one can see, this is the same Generalized DMPK equa- 
tion as used in Ref. Q- But this time we do not need 
to assume anything special about the eigenvector corre- 
lations, and indeed by including them, the equation now 
automatically conserves probability. Finally, since this 
equation was derived under very general assumptions, we 
expect it to provide a description of the conductance dis- 
tribution on both sides of and across the metal-insulator 
transition. The only remaining problem is that Eq. (j2 . 14[) 
implies that the level repulsion for neighboring levels is 
determined by the matrix 7 {, which is not obviously re- 
lated to the symmetry parameter ft. In the insulating 
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regime, it is known that on average jab <C 1 which seems 
to contradict this symmetry requirement. In the next 
section we address this issue. We show that under reason- 
able assumptions, 7 a fc — > ft in the small separation limit. 
To demonstrate this we develop an evolution equation for 
the matrix K a b- 



mean value of k in the weakly disordered quasi-ld regime, 
and also in the 3d metallic regime, as shown in Ref. [2j. 
Although these fluctuations are not generally small in 
the 3d strongly disordered regime, we may expect them 
to be negligible when the eigenvalues assume metallic-like 
configurations, namely, 1 — A„/A n +i <C 1/N. For such 
levels, Eq. (|3.6I) will reduce to: 



III. LEVEL REPULSION BETWEEN NEAREST 
NEIGHBORS 

We begin developing an evolution equation for K mn by 
expanding K mn out to linear order in 8L z /l. We write 



■5K, 



<nn - £ (\ Vma + Sv ma\ \v na + 5v na | 2 ) (3. 1) 



where the brackets denote an average over v. Then using 
the results of Eq. (12.91) . Eq. (|2.11[) . and the random phase 
assumption we find that we can write the evolution of 
K m n with length L z as 



dK„ 



dt 



E f ( j mi j irn \ , \ y r ( j ni j i 



i^m 



i^n 



r\ \ ^ r j mi £ r\ r j rim c 

* / Jim- 1 -' mi ran ^Jmn^-'nm m^n 

i^m 



where /y is given by, 



_ AiCl + AO + Ajd + Aj) 



(Aj - X 3 ) 2 



(3.2) 



(3.3) 



and L J ik is given by, 



Lf k = ^(\v ia \ 2 \v ka \ 2 \v jf3 \ 2 \v kfi \ 2 ) 

a,f) 

= (hkkjk)- (3.4) 

The evolution equation for K conserves probability for if 
we start with a K mn whose columns sum to 1, as it must, 
then one can show that the evolution equation preserves 
this property. One can also see from equation Eq. (I3.2|) 
that the evolution of K a i, = {k ao ) is determined purely 
by its second moment. To gain some further insight we 
will split L into the sum of two terms, its mean field 
approximation and fluctuations about the mean: 



L 



jk _ 
ik 



(k l k)(k 3 k)+&H=K lk K ] k+^. (3.5) 



Next we set m = n in Eq. (13.2[) . and combine it with 
Eq.daH]) to obtain: 



dK 



dt 



mm _ jr2 

mm 



7mi(l - 7„ u ) + 



2A™ 

mi 



A Arm 
mm 



K 2 

mm 



This equation has a few interesting consequences. First, 
when fluctuations in k are small, we may neglect the 
A term. Such fluctuations will be small compared to the 



dKr, 



,2 \ " 

v mm / j . 

iy^m 



im ^imi ( 1 *Ymi ) • 



(3.7) 



Now if we let nearest neighbors A m and X n approach each 
other in Eq. p.7[) , f mn will become singular. But since 
dKmm/dt must be bounded, it must be that 7 mrl goes 
to I. And so we see that in the small separation limit 
|A„ — A m | — > 0, 7„ m ft = 1, as required. Thus for 
instance, even when Ai and A2 are far apart on average 
and consequently 712 < 1 in the large disorder regime 
(as shown in numerical studies 2 ), in the limit where their 
separation goes to zero, the repulsion will be determined 
by the usual symmetry parameter equal to ft = 1 in the 
orthogonal case as considered here. 

Unfortunately it is difficult to access the above limit 
numerically since for large disorder the lowest levels 
have a statistically negligible probability to become close. 
However it turns out that there are some levels near the 
middle of the spectrum that do come closer for any given 
disorder and for sufficiently large system size (see Fig. 
ID- Fig. [T] illustrates the trend of how 7, ljn +i changes 
when the separation of the levels becomes smaller. Here 
we have used the same anisotropic Anderson model as 
considered in Ref. • We consider a cubic system of size 
Lx Lx L but with anisotropic nearest neighbor hopping 
elements t z = 1 and t x = t y = 0.4. The anisotropy guar- 
antees that all channels are open, i.e. N — L 2 . For a 
given disorder W and a fixed system size L, the matrix 
K and hence the parameters 7 mn as well as the eigenval- 
ues A n are then evaluated using 10 4 statistical samples 
within orthogonal symmetry. 

Since the eignevalues A„ = cosh _2 (a; n /2) decrease ex- 
ponentially with the size of the system, it turns out that 
many eigenvalues decrease below our numerical accuracy. 
Therefore, we consider only channels where (x n — xi) < 
34. This guarantees that A„/Ai > 1.7 x 10~ 15 , although 
at the same time the restriction does not enable us to 
consider all 7 mn for a given N. For instance, in Fig. [TJ 
n < 120. Nevertheless, Fig. [T] confirms our expectation 
that "fnm+i increases towards unity when the difference 
(x n+ i — x n ) decreases. 

On the other hand, we expect Eq. (|3.7p to also describe 
a weakly disordered quasi-ld conductor. In the large 
L z limit, the eigenvector correlations approach a limiting 
distribution. In this case we must have dKmm/dt = 0. 
And it follows that 7 mn = 1 as well, in accordance with 
the DMPK equation. Indeed, if fluctuations are negli- 
gible we can solve for K directly from Eq. (|3.2j) . Using 
Eq. (I3.5[) with A = 0, the evolution equation for K re- 
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Xin+1 



0.1 



W= 16 




o.i - 



4 % -/V 




j i i |_ 



0.1 



< x , - x > 

n+ 1 n 



FIG. 1: Color online: 7n,n+i as a function of the difference 
(x n +i—x n ) for statistical ensembles of cubic samples of differ- 
ent sizes 10 < L < 30 for a given disorder strength W = 16. 
(Only N = 100 samples are considered for L = 30.) For com- 
parison, the critical disordes^ W c ~ 9.4. To avoid numerical 
inaccuracies, only channels with {x n — xi) < 34 were consid- 
ered. The mean values of the logarithm of the conductance 
are -4.07 for L = 10 and -12.6 for L = 30. Typically, the 
level separation decreases with increasing n up to n = N/2, 
and then start to increase again; this leads to two different 
'branches' (the upper one for n < N/2), as is visible for the 
two smallest sizes L < 14. For larger L this crossover occurs 
at larger n and smaller separation, the lower branch eventu- 
ally becoming numerically inaccessible for sufficiently large L. 
Inset: The same ( r f n ,n+l vs the level separation) for a fixed 
length L — 22 but for different values of W. 



duces to: 



dt 



+ 2) fimK^iS, 



K 



nm- 1 - l tm) 



mn-** nra th^7i 



(3.8) 



It is straightforward to show that the L z independent 
solution to this equation is 



K — 



1 ~T" &mn 

N + l 



(3.9) 



as required by the quasi Id DMPK equation. In order 
to move beyond the small fluctuation limit, and model 
more generally a strongly disordered insulator, especially 
near the critical point, it is expected that A will become 
important. Though little can be said at this point, one 



implication of Eq. (|3.6[) is that if K rnm reaches a limit- 
ing distribution in L z L strongly disordered limit, the 
steady state value of -f mn will depend on these fluctua- 
tions. 



IV. THE K-MATRIX AND THE LYAPUNOV 
EXPONENTS 

Finally, we develop an equation relating the Lyapunov 
exponents and the matrix K. The Lyapunov exponent v 
is defined by 

\(L z ^f oo) = e 2vL *. (4.1) 
In order to develop an equation for v we write 



X(L Z + SL Z ) - \{L Z ) = 
Solving for v we have 



_ 2v(L z +5L z 



„2vL, 



25L, n V + A 



(4.2) 



(4.3) 



Expanding out to linear order in SL z /l and averaging 
over disorder we obtain 



2SL Z 




Using Eqs. (|2.9p and (|2.11l) . we finally arrive at 
1 



-K n 



V 1 + Xn K 



An — An 



(4.4) 



(4.5) 



for all disorder. A similar equation was obtained in 
Ref. |l9j . Eq. f)4.5[) has been verified numerically in the 
insulating regime, where most doubt concerning its va- 
lidity would reside. When plots of v n obtained via nu- 
merical simulations using tight binding Anderson model 
and from Eq. (I4.5[) are overlaid, as shown in Figure [21 the 
points are indistinguishable. 

Interestingly, if we sum this equation over n, we obtain 
the following relationship, independent of K and A, 



" N 



(4.6) 



This relation has been known to be true in the diffusive 
limits, but the present calculation implies it is true for 
large disorder as well. This has been confirmed numeri- 
cally. 

Note that these relations are independent of any model 
for the matrix K . A 2-parameter model for K was chosen 
in Ref. [H for large disorder, based entirely on numeri- 
cal studies. The exact relationship should allow us to 
build a better model based on the known properties of 
the Lyapunov exponents. 
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FIG. 2: (Color online) The Lyapunov exponents v n plotted 
vs. n. The red diamonds are from numerical simulations 
(L = 10, W = 20, and <ln g> = -7.1); the solid blue circles 
are from Eq. (I4.5|l . 



V. SUMMARY AND CONCLUSION 

In this work we developed a systematic and controlled 
derivation of the Generalized DMPK equation by ex- 
panding the probability distribution Pl^+slA^) m a 
power series 8L z /l and keeping all terms linear in SL z /l 
in the limit SL z /l — > 0. This means including certain 
eigenvector correlations that had been neglected before. 
The additional correlations do not change the Gener- 



alized DMPK equation obtained earlier; instead, their 
inclusion automatically conserve probability by allowing 
the evolution of the Jacobian without having to rede- 
fine it with a disorder dependent exponent. In addition 
to providing a broader applicability of the Generalized 
DMPK equation such as the conductance distributions 
obtained in Ref. [l|, the derivation shows how the re- 
pulsion of closely spaced neighboring eigenvalues remain 
consistent with symmetry requirements, even though it 
can be very different when the eigenvalues are far apart. 
Moreover, we obtain an exact result relating the correla- 
tion matrix K and the eigenvalues A with the Lyapunov 
exponents v as well as a sum rule for the exponents in- 
dependent of K or A. 

The challenge remains to construct a reasonable model 
for the phcnomcnological matrix K. The evolution of K 
involves quantities that can not be factored into products 
of -ftT's in general, and at this point an analytic solution 
seems too complicated. On the other hand a crude phe- 
nomenological 2-parameter model suggested by numer- 
ical studies and used in Refs. [lH seems to work quite 
well in the strongly disordered regime, although further 
numerical studies are needed to check if /how the model 
might change when e.g. approaching the metal-insulator 
critical point. We hope that the relationships with the 
Lyapunov exponents obtained here should provide useful 
insights. 

PM acknowledges financial support from the Slovak 
Research and Development Agency under the contract 
No. APVV-0108-11. 
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